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We  rJescribe  the  change.s  in  refractive  index  induced  i>y  the  photoelantic 
effect  when  an  optical  fiber  in  subjected  to  a uniformly  applied  diametral 
ntressy  F’or  easily  achievable  values  of  the  force  per  unit  length  applied  to  the 
fiber^+Ciue  find  that  regions  of  equal  or  higher  refractive  index  than  the  core 
may  be  induced  in  the  outer  region  of  the  fiber.  Thun  the  ntrensed  region  is 
capable  of  acting  as  a mode  converter  that  affects  the  transmission  character- 
istics of  the  fiber. 
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Refractive  Index  Changes  in 
Optical  Fibers  Subject  to  Diametral  Stress 


I.  iMiionn  TioN 

When  an  optical  fiber  is  subjected  to  a uniform  stress  applied  along  its  diam- 
eter, refractive  index  changes  will  be  induced  owing  to  photoelasticity.  This 
phenomenon,  one  of  several  which  can  alter  the  transmission  through  optical  fibers, 
may  be  involved  in  the  observed  changes  of  the  transmitted  mode  patterns  in  exper- 
iments conducted  at  Solid  State  Sciences  Oivisi<m  (h'TSS).^  Assuming  that  a uniform 
stress  is  applied  to  a perfect  fiber,  one  may  also  induce  changes  in  the  geometri- 
cal cross  section.  In  the  nonuniform  case,  which  may  be  closer  to  reality  in  prac- 
tice, one  may  induce  microbending  and  inhomogeneously  stressed  regions,  both  of 
which  are  capable  of  altering  optical  transmission.  The  present  work  is  limited 
to  an  examination  of  the  changes  in  refractive  index  and,  hence,  the  waveguiding 
characteristics  of  a fiber,  induced  by  the  photoelastic  effect. 


(Heceived  for  publication  1!)  April  1977) 

The  term  uniform  is  used  in  the  sense  of  how  the  external  force  is  applied  and  the 
geometrical  shape  of  the  fiber, 

^The  authors  thank  Dr.  H,  Payne  for  discussit>ns  of  his  experimental  results. 
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2.  sti«kss-im)h:i:ii  ( II\\(.ks  in  ki;i  H\r;n\i;  imma 


Whenever  an  optical  medium  id  stredsed,  itd  refractive  index  becomed  altered, 
a phenomenon  which  id  termed  the  photoeladtic  effect.  If  we  represent  the  refrac- 
tive index  by  n,  the  streds  by  o,  and  the  I’ockeld  piezooptic  constants  by  q,  then 
the  change  in  reciprocal  of  the  refractive  index  squared  is  given  as:^ 

^ a (k,  j 1,  fJ)  , (1) 

n'^  J J 

o 


where  the  standard  (Voigt)  compact -subscript  notation  is  used,  and  n^  is  the 
refractive  index  of  the  unperturbed  medium. 

Let  us  set  up  a coordinate  system  such  that  the  z-axis  is  parallel  to  the  cylin- 
drical guide  axis,  the  y-axis  parallel  to  the  direction  along  which  opposing  external 
diametral  forces  are  applied,  and  the  origin  x y 0 is  located  at  the  center  of  the 
guide.  The  forces  (F)  extend  for  a specified  length  (I.)  along  the  z-axis  (see  I’ig- 
ure  la).  With  this  setup,  any  arbitrary  point  (x,  y)  in  the  fiber  experiences  tensile 
or  compressive  Stresses  parallel  to  the  x-  and  y-axes  (represented  by  Oj  and 

respectively,  and  shown  in  Figure  lb)  and  shear  stresses  o {=  a ) which  tend  to 

t)  xy 

cause  rotation  about  the  z-axis.  However,  this  latter  stress  vanishes  on  the  x- 
and  y-axes.  Under  these  circumstances,  the  total  fractional  change  in  refractive 
index  due  to  the  stress  at  the  point  (x,  y)  is  given  by: 


2 

-n 

(qil  - - <^2*x,  y 


(2) 


for  an  optically  isotropic  material.  We  initially  assume  that  the  force  F is  applied 
at  the  fiber's  surface  by  knife  edges  running  parallel  to  the  z-axis  for  a length  L. 

The  general  expressions  for  the  quantities  y and  (o2^x  y arbitrary 

point  (x,  y)  are  derived  by  Frocht.  ^ Along  the  x-axis  they  are  given  by: 


«^l>x.0  = » 


^2 

2 

D 

- 4 

X 

2 

2 

D 

t-  4 

X 

(3) 


1.  Nye,  J.  F.  (1957)  Physical  Properties  of  Crystals,  Oxford  Press,  London, 

p.  224. 

2.  Frocht,  M.  (1948)  Photoelasticity,  Vol.  2,  Wiley,  New  York,  pp.  125-129. 

II 


(0  ) (b  ) 


I-'igure  1.  Coordinate  System  for  the  I'iber,  Showing  tiie  Applied 
External  Eorces  (!•)  and  the  Hesnlting  Internal  Stresses  (o)  at  an 
Arbitrary  Point 


(ct„)  _ 

2 X.  0 


-II 


4 D 


4 X ) 


- 1 


where  D is  the  fiber  diameter,  the  range  of  x is  0 ■?  ] x|  5"  | D/2I , and,  the  quantity 
H is  defined  as; 


\ tD  I L • 


Along  the  y-axis  (0  ^ | y|  | 0/2]  ),  these  two  stresses  become: 


y 


-HD 


n - 2 y D ( 2 y D 


(o) 


Hi) 


(7) 


The  negative  signs  in  the  a's  indicate  that  they  are  compressive  stress.  The 

apimoximato  behavior  of  Eqs.  (3),  (4),  (h)  and  (7)  is  shown  in  h'igure  2.  These 

curves  are,  of  course,  sytnn^etrical  about  the  origin.  On  each  plot  the  positive 

vertical  axes  pertain  to  n whereas  the  negative  vertical  axes  refer  to  a^.  Note 

that  the  (o,,)„  curve  in  I'igurc  2b  becomes  infinite  at  the  endpoints  (y  tD/2). 

" tt,  y 

This  occurs  because  the  force  is  being  applied  over  an  infinitesimally  small  surface 


area. 
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1 a ) 


( b) 


2.  Stres.ses  in  I*’iber.  and  02  distance  along  the 
horizontal  (x)  axis  and  vertical  (y)  axis  measured  from  the  fiber's 
center 


When  the  foregoing  information  is  substituted  into  Eq.  (2),  we  g'  after  inter- 
changing with  qj  j : 


(“) 


■S  n^ 

(qi2  - ■ °2^x,  0 " “ *‘^12  ■ '’ll' 


X.  0 


(4  I1U^(D^  - 4 x^)/(D^  ♦ 4 x^)^)  , 


(8) 


and 


2 2 

-"aVy  ^lq.2-<l„>|4HD'';(D2-4/)l 


(9) 


At  the  center  of  the  cylinder  (x  y 0),  both  of  the  preceding  equations  reduce  to; 
2 _2 


(•^)p  0 'T' *'’12  ■ ‘’ll'^'’l  ■ '’2'o.  0 2“ ’'>12  ■ '’ll* 


"o  *‘’12  ■'’11*  *7tD  ■ L 


(10) 


8 


Kquation  (8)  shows  that  (^n  n)^  has  its  maximum  value  at  x o,  rliminishe.T 
witli  X m almost  the  same  manner  as  Kq.  (8)  (see  also  h'igure  2a),  then  vanislies 
at  X il)  2.  On  the  other  hand,  Eq.  (9)  behaves  (|uite  similarly  to  Eq.  <7)  for 
eliancing  y,  including  the  singularity  at  the  endpoints. 

However,  if  E is  applied  by  means  of  flat  (dates,  I'ather  tlian  knife-edges,  so 
that  the  force  is  distidhuted  over  more  of  the  filler's  surface,  then  the  (loints  of 
maximum  stress  move  a short  distance  away  from  the  edges,  closei-  to  the  fiber's 
center.'  We  designate  the  positions  of  these  [loints  by  ty^j.  Since  the  force-per- 
unit  length  enqiloyed  in  our  case  is  of  the  same  order  of  magnitude  as  in  the  exam- 
ple treated  by  I'rocht.'^  we  should  expect  similar  qualitative  results,  namely,  that 
the  maximum  stresses  along  the  y-axis  move  to  a position  ap()roximately  tD  'sO 
away  from  the  edge  and  that  the  quantity  (a,  - cj,p)„  be  a[)[)roximately  ten  times 
greater  than  that  at  the  center,  that  is. 


°2'o,  y ' ■ <^2^0,  0 


As  a consequence. 


At  the  edges  (y  tD/2),  the  magnitude  of  (a^  - ,,,  and  (.in  ii)^  ,,  then  fall 

to  approximately  three  times  the  corresponding  values  at  the  center  of  the  fiber. 

The  term  (ntj  - ^ is  plotted  in  Figure  ,'la  as  a function  of  x.  It  is  given  by 

the  quantity  inside  the  curly  braces  on  the  right-hand  side  of  Ec(.  (8).  As  stated 
previously,  it  behaves  somewhat  like  the  (Cj)^  q curve  of  f''igur('  2a.  In  I'igure  8b, 
however,  we  have  plotted  two  versions  of  the  quantity  'I'lie  solid 

curve  pertains  to  the  case  in  which  the  forces  are  applied  by  knife-edges.  It  is 
governed  by  the  term  inside  the  curly  braces  on  the  right-hand  side  of  Eq.  (9). 

The  dashed  curve  pertains  to  the  case  in  which  the  forces  are  ap|ilied  by  flat  (dates. 
In  accordance  with  Frocht,  the  profiles  for  both  cases  have  the  same  magnitude, 
namely,  4 II,  at  the  center  and  behave  identically  out  to  a distance  of  about  n 10 
from  the  edges.  Thereupon,  the  dashed  curve  peaks  at  a value  ' 10  times  that  at 
the  confer  (^  40  II)  at  a point  ' D,  50  from  the  edges  and  then  falls  to  ' 12  II  at  the 
cfiges. 


8.  See  Kef.  2,  [>p.  27-29. 


Figure  Stress  Differences  in  ('"iber.  (a^  -02) 
vs  distance  along  x-  and  y-axes.  'i'fie  solifi 
curves  pertain  to  application  of  forces  via  knife- 
edges;  the  dashed  curve  in  (b)  via  flat  plates. 


We  now  fxirticularize  our  treatment  to  the  case  of  a fiber  having  an  outer 
diameter  (D)  of  5 mils  and  a core  diameter  of  3 mils.  [ F’or  this  fiber  we  get 
11  137  b'  I . (Ib  'ii/)  from  Fq.  (f)).  | The  fiber  is  assumed  to  have  a stepped  index 

profile  in  the  unstressed  state,  the  difference  in  indices  between  core  and  cladding 
being  1 percent.  Designating  this  unstressed  refractive  index  change  by  (.An  n)  , 
we  plot  this  versus  distance  along  any  arbitrary  diameter  in  Figure  4a. 


If  we  assume  that  the  fiber  core  is  fused  silica,  then  n 

1-S  - 


1.  a and  - q^  ^ 


2 10  in  /lb  (see  Hef.  4).  If  we  consider  the  case  in  which  the  forces  are 

applied  by  means  of  flat  plates,  then  we  could  calculate  the  specific  values  of 
(An/n)''’  along  the  x-axis  via  Eq.  (8)  and  along  the  y-axis,  using  Fqs.  (10)  and  (12), 
filling  in  the  rest  of  the  values  by  utilizing  the  information  given  previously. 

(See,  for  example,  the  information  contained  in  I'igure  3b.)  Tlie  total  i-efracti\’e 

T 

index  change,  designated  by  (.An/n)  , would  then  he  given  by  the  sum: 

y 


(An/n) 


X,  y 


(An/n)|^  y ' (An/n)^  ^ 


(13) 


4.  Kapron,  F.  ( 1 952)  Rirefringence  in  dielectric  optical  waveguides,  IKKJO.l.  of 
Q.  Elect.  QE-«:224. 
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( b ) 


' 0 I m I I 


(AJL 
V n 


I'igure  4.  ({efractive  Index  Change  v.s  Distance:  (a)  unstressed 
ease;  (b)  stress  alone;  and  (c)  total  effect 

SI' 

In  Table  1 we  list  the  (An/n)'  and  (An  n)  for  the  above  silica  fiVjer  at  3 loca- 
tions, namely,  the  center  and  the  two  symmetrical  (joints  (0.  ty,,)  for  various 
values  of  T 1,.  (.Mote  that  (An /n)  is  alway.s  0.  01  at  the  center  and  zero  at  the 
points  (0,  iryj^j))  . 

In  I-igures  4b  and  c we  utilize  some  of  the  information  contained  in  Table  1. 

In  these  two  figures,  v.e  plot  the  index  changes  due  to  stress  alone  and  the  total 
effect,  restJectively,  along  the  y-axis  for  F/I,  - 100  and  200  Ib/in.  .Some  of  the 
important  features  of  the  curves  in  Figure  4c  are  brought  out  much  more  vividly 
in  f igure  .'j.  Here,  we  have  plotted  the  total  index  change  as  a function  of  F/L  for 
3 positions  in  the  fiber,  namely,  (0,0)  and  (0,  Note  how  the  total  index 

change  at  the  center  barely  increases  in  magnitude  with  change  of  applied  F/I,. 

On  the  other  hand,  the  total  index  change  at  the  point  of  maximum  stress  turns 
out  to  be  a very  sensitive  function  of  F/L.  At  F/L  50  Ib/in,  its  magnitude  is 
slightly  greater  than  half  that  at  the  fiber's  center;  at  100  Ib/in.  it  is  approximately 
the  same;  while  at  200  Ib/in.,  it  becomes  almost  twice  that  at  the  center. 

1 1 
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'J  able  1.  Kefraetive-  Index  Chanfje  vs.  I'<jrce  Applied 


It  .-.lioiild  he  l ememberefl  tliat  tlie  re.-iults  di -iplayed  in  I igures  4h,  4c-  and  5 
hold  only  along  the  y-axi.s.  As  c-an  he  seen  with  the  hc-lp  of  Kcp  (ti)  |and  1 igure  .'la| 
the  (^n  n)^  profile  along  the  x-axis  sIkjws  a maxirmim  at  (0,  0)  and  clrops  off  to 
zero  at  the  [leriphery  (x  tIV2).  However,  even  this  niaxinium  is  still  one  order 
c.r  magnitude  lower  than  the  original  value  of  (An,n)^  ' O.Ol.  Consecnient ly,  the 
ront  rihution  to  An  n along  the  x-axis  due  to  th<‘  stit-ss  is  unnnportant  and  its  total 
index  profile  would  still  he  given  as  shown  in  l igure  4a.  I-'or  diametral  lines  at 
arhifrary  angles  with  respect  to  the  x-  and  y-axes,  the  (An  n)  * profiles  would  he 
somewhere  between  those  shown  in  1-igures  4a  and  e. 

When  the  applied  linear  forces  are  such  as  to  cause  large  index  changes,  as 
shown  m Figure  4c,  for  example,  then  one  must  take  into  consideration  the  possi- 
bility of  coupling  at  least  some  of  the  optical  signal  from  the  core  to  *he  high  index 
region  near  the  periphery.  The  net  effect  of  this  would  he  the  conversion  of  modes 
in  the  core  to  the  stressed  region,  coupling  of  energy  out  of  the  guide,  and  deteri- 
oration of  the  original  signal.  This  process  of  optical  coupling  would  ho  analogous 
to  quantum  mechanical  tunneling  between  two  potential  wells  separated  by  a barrier. 

As  a convenience  for  the  reader,  we  have  included  two  appendices.  In 
Appendix  A w-e  list  the  material  parameters  for  fused  silica  which  would  be  appro- 
priate for  I notoelastic  analyses.  In  Appendix  B we  have  calculated  the  strains  at 
the  points  (0,  0),  (0,  <0  tD/2)  and  (tD/2,0)  for  the  particular  fiber  treated 

in  this  paper.  , 


:l.  ( (INU.I  klONS 

We  have  shown  that  a fiber  subjected  to  easily  achievable  values  of  applied 
stress  will  assume  a refractive  index  distribution  in  which  portions  of  the  outer 
region  (usually  the  cladding)  possess  an  equal  or  greater  refractive  index  than 
that  of  the  active  waveguiding  region  (core).  Hence,  the  stressed  region  is  capable 
of  supporting  new  modes  and  thereby  altering  the  distribution  of  the  old  ones;  that 
is,  it  is  capable  of  acting  as  a mode  mixer  or  converter.  The  internal  stresses 
generated  by  the  applied  force  are  typically  an  order  of  magnitude  or  more  smaller 
than  those  corresponding  to  the  intrinsic  strength  of  fused  silica.  Although  a 
detailed  calculation  of  the  mode  conversion  induced  by  transmission  through  the 
stressed  region  would  be  prohibitively  difficult,  it  nevertheless  appears  safe  to 
conclude  that  the  mechanism  described  here  can  lead  to  perceptible  changes  in 
fiber  transmission. 
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Appendix  A 

Fused  Silica  Parometers  Pertinent  to  Photoelastieif)r 


The  following  quantities  are  taken  from  the  AlP  Handbook^  and  converted  into 
English  units: 

E = Young's  Modulus  = 1.06  y 10^  Ib/in^ 
r = Poisson  ratio  = 0.  17 

7 2 

H - Larne's  Constant  - .453  y 10  lb/ in 

7 2 

X = Lame's  Constant  - .234  y 10  Ib/in 

2 

The  values  of  Pockels  elastooptic  roefficients  are  taken  from  Pinnow 


11 


121 


Pl2  = .270 

The  following  values  of  elastic  compliances  stiffness  constants  (c^^)  and 

F^ockels  piezooptic  coefficients  (q^^)  are  computed  from  the  above  quantities: 


1.  Gray,  D.  E. , Ed.  (1972)  AIP  Handbook,  3rded.,  McGraw-Hill,  New  York, 

p 3-104. 

2.  Pinnow,  D.A.  (1972)  Electro-optic  materials,  in  Laser  Handbook.  Voi.  1, 

F.  Arecchi  and  E.  Schulz-Dubois,  Ed.,  North  Holland  Publ.  Co., 
Amsterdam,  p.  999. 
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Appendix  B 

Calculation  of  Strains 


Strains  are  related  to  stresses  (a.)  in  a material  medium  via  the  elastic 
compliances  (s.^)  in  the  followinfi  manner: 

f-  s.  o.  (i,j  1,  fi)  . 

> 1.1  J 

Again,  the  compact  subscript  notation  has  been  employed.  However,  along  the 
X-  and  y-axes,  only  Oj  and  exist.  Hence,  the  x-  and  y-  components  of  strain 
along  these  axes  in  an  isotropic  material  can  be  written  as: 


*^1  ’ 1 1 ” 1 1 2 °2 


* 2 '*12  ”l  ■ ■‘^1  1 '’2  • 


t sing  the  same  fused  silica  fiber  as  mentioned  in  the  text,  in  which  n = 1.  5, 

-7  2 " 

I)  .00.')  in,  S|j  and  Sj^  are  0.  9S  and  -0.  10  (x  10  in  /lb),  respectively,  with 
I-  1,  100  It) 'in,  we  calculated  the  following: 


At  X y 0: 


1. 8 

-3.  8 


X 10 
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X 10 
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X 10'^ 

because  = 02  - 0 


18 


METRIC  SYSTEM 


BASE  UNITS 

Quantity 


Unit 


IcnKth 

ma.s« 

tim«* 

eln  tri(  c urrpnl 
themiixlvnamu  temperature 
amount  of  substani  e 
luminous  intensity 

Sl'PPLEMENTARY  I'NITS 

plane  an^^le 
solid  angle 

DFJIIVIJ)  I NITS 
A<  I eleration 

a<  tivitv  (of  a radioai  tive  soun  e) 

angular  art  eleration 

angular  velocity 

area 

density 

elet  tri«  (.aparitanie 

elet  tri<  al  rnndu<  lani  e 

elet  trie  field  strength 

elettrii  inductance 

elet  tnt  potential  difference 

elet  trii  resistant  e 

elet  Iromtitive  fort  e 

energy 

entropy 

fort  e 

fretjuem  y 

illuminant  e 

luminance 

luminous  flux 

magnetit  field  strength 

magnetK  flux 

magneiK  flux  density 

magnetomtitive  forte 

power 

pressure 

(quantity  of  elet  tnt  ity 
f)uanliiv  of  heat 
radiant  intensity 
spet  ifir  heat 
stress 

thermal  ttmdut  livity 
velfK  itv 

visfoaity  dynamo 

vtst  iMity  kinematic 

vtiltage 

volume 

wavenumber 

work 


metre 

kilogram 

set  ond 

ampere 

kelvin 

mtile 

( andela 


radian 

sleradian 


metre  per  s«.ond  stjuared 

disintegratitin  per  set  ond 

radian  per  set.ond  squaretl 

radian  per  settmd 

square  metre 

kiltigram  per  t.ubit  metre 

farad 

siemens 

volt  per  metre 

henry 

volt 

ohm 

volt 

K>ule 

fouJe  per  kelvin 

newton 

hertz 

lux 

t.andela  per  square  metre 
lumen 

ampere  per  metre 

weber 

tesla 

ampere 

watt 

(Mst  al 

( oulomb 

|oule 

watt  per  sleradian 
joule  per  kilogram-kelvin 
past  ai 

watt  per  metre  kelvin 
metre  per  set  onti 
past  al-set  tmd 
s(|uare  metre  per  set  ond 
volt 

cubit  metre 
renprr>r.al  metre 
joule 


SI  PREFIXES 


Multiplication  Factors 

1 non  (KMMKKMMH)  10'^ 

1 OOtMNHMKNl  1(1* 

I IKMI  000  • 10* 

I 000  « 10* 
100  - UH 
1(1  • 10* 

0 1 1(1  ' 

0 01  10  * 
0 001  • to  * 

OfHIOOOl  10  * 
0 (KNI  (NHI  001  " H)  * 
0 000  (MK)  0(M>  001  « 1(1 
O 000  000  (KNI  OOO  001  10 

(XHKMNlOOOOtNKXNNIOOl  10  '• 


SI  Symbtd  Formula 

m 

kg 

s 

A 

K 

mol 
< d 


rad 

$r 


I 

S 

H 

V 

V 

) 

N 

Hz 

lx 

Im 

VVb 

T 

A 

V\ 

Pa 

C 

1 


ms 

(disintpgralionj's 
rad  s 
rads 
m 

kgm 

As-V 

A\’ 

\ m 
Vs.  A 
V\  A 
V A 
W A 
N*m 
I'K 

kgms 

(cyi  le)/s 

Im'm 

edm 

cd*sr 

Am 

Vs 

Wbcm 


)'» 

N m 
A*s 
N-m 
VV  sr 
Ikg-K 
N m 
W mk 
m s 
Pas 
m s 
W A 
m 

(wave)m 

N-m 


I'rnfu  M Symt«'l 


lera 

K>H* 

mega 

kilo 
hiH  to* 

deka* 

det.l* 
t entl* 
mllli 
mil  n> 
nano 
pit  II 
femto 
attic 


r 

(; 

M 

k 

h 

da 

d 

I 

Ml 

n 
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MISSION 

of 

Rome  Air  Development  Center 


RADC  plans  and  conducts  research,  exploratory  and  advanced 
development  programs  in  command,  control,  and  conmunications 
(C^)  activities,  and  in  the  areas  of  information  sciences 
and  intelligence.  The  principal  technical  mission  areas 
are  communications,  electromagnetic  guidance  and  control, 
surveillance  of  ground  and  aerospace  objects,  intelligence 
data  collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microvave 
physics  and  electronic  reliability,  maintainability  and 
compatibility . 


Printed  by 

United  State*  Air  Force 
Honecom  AFB,  Mo**.  01731 


